UFIFT-HEP-00-17 



hep-ph/0007019 



Analyzing the 't Hooft Model on an x^-p'^ Lattice^ 



Joel S. RozowskyQ and Charles B. ThornQ 

Institute for Fundamental Theory 
Department of Physics, University of Florida, Gainesville, FL 32611 

(February 1, 2008) 



Abstract 

We study the 't Hooft model (large Nc QCD in 2 space-time dimensions) using an improved 
approach to digitizing the sum of gauge theory Feynman diagrams based on light-cone gauge 
— and discretized and ix^ . Our purpose is to test the new formalism in a solvable 
case, with the hope to gain some insight into how it might be usefully applied to the physically 
interesting case of 4 dimensional QCD. 
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1 Introduction 



Last year, with Bering we proposed fll a new method to digitize the sum of planar diagrams selected 
by 't Hooft's Nc ^ oo limit of SU{Nc) gauge theories 0]. The proposal, based on the light-cone 
or infinite momentum frame description of the dynamics, involved discretization of both the 
carried by each line of the diagram and the propagation time r = ix~^ , as in 1^-^ . But the main 
advantage of the new version was a coherent prescription for resolving most of the ambiguities due 
to p+ = divergences that typically plague the light-cone description. 

We hope that our formalism will eventually allow an improved understanding of QCD in 4 
dimensional space-time. But in this article, we merely wish to test the proposal in the context of 
the well- understood case of large Nc gauge theories in two space-time dimensions, namely the 't 
Hooft model p. Our purpose is not to unearth new aspects of the model, but rather to see how 
its well known properties can be obtained from our new discretization. 

The physical content of the 't Hooft model boils down to an integral equation, essentially a 
Bethe-Salpeter equation Q, that determines the mass spectrum of qq mesons. The reason the 
limit A'c ^ oo reduces to ladder diagrams (albeit with self-energy corrected quark propagators), 
is that the 2 dimensional gluon is not dynamical (there are no transverse polarizations). Thus, 
as with any axial gauge, the light-cone gauge A- = eliminates all gluon self-interactions, so 
j4_|_ can be integrated out inducing an instantaneous Coulomb potential. But the 't Hooft limit 
Nc ^ oo further eliminates all quark loops and all non-planar diagrams, leaving only the planar 
self energy corrections to the quark propagator, and the ladder bare gluon exchanges (Coulomb 
interaction) between quark anti-quark lines in the singlet qq channel. In light-cone parameters the 
Bethe-Salpeter equation summing these ladder qq diagrams simplifies to the single variable 't Hooft 
integral equation Q. 



M 'Pi^)={z + Y^)f^^i-)-^Pl dy , (1.1) 



where the integral is understood to be evaluated by the principal value prescription. The variable x 
is the fraction carried by the quark of the total -P"*" of the system (the anti-quark carries fraction 
1 — x). Also A4 is the mass of the meson bound state and if satisfies the boundary conditions, 
ifiO) = ifil) = 0. 

Since the new formalism discretizes r = ix'^ = ka in addition to p"*" = Im, the corresponding 
simplifications lead to an equation that is not a straightforward discretization of this integral 
equation. In particular, the continuum limit can be taken in different ways depending on the ratio 
Tq = m/a (which would be infinite for continuous r), and we want to explore to what extent these 
different continuum limits lead to the same physics. We shall find that some care must be taken 
with the setup of the discrete r dynamics in order for this to be true. Indeed, a numerical study 
shows that the most simple-minded treatment leads to a ground state that becomes unstable at 
moderate 't Hooft coupling even with relatively small P~^/m = M unless the ratio a/m = I/Tq is 
tuned to be sufficiently small (perhaps infinitesimal for large M). If this feature were robust, it 
would cast doubt on any potential utility of the discretization of r. 

To overcome this difficulty, we find it necessary to veto some of the "densest" discretized Feyn- 
man diagrams: a quark must be forbidden to emit 2 gluons at immediately successive time steps, 
with a similar veto on two successive absorptions. With this simple veto (which is prescribed locally 
in time), we shall show that the continuum limit reduces to the 't Hooft model provided only that 
the total P^ of the qq system is large compared to the discretization unit m. In particular it is 
not necessary that the ratio Tq = m/a be large. Keeping Tq finite in the continuum limit leads to 



1 



the 't Hooft equation with a non-trivial renormahzation of the couphng. Because of this effect, it 
turns out that the effective (renormahzed) couphng is small for both large and small bare coupling, 
reminiscent of strong/weak coupling duality. The strong coupling limit favors the densest diagrams, 
so vetoing some of the densest ones has a dramatic effect on the strong coupling behavior of the 
theory. This possibility was anticipated and discussed in Q in connection with the nature of the 
fishnet diagrams in higher dimensional space-time. 

The rest of the paper is organized as follows. In Section 2 we set up the discretized 't Hooft 
model. We analyze it using a single time-step transfer matrix in Section 3 and using a Bethe- 
Salpeter approach in Section 4. In Section 5 we discuss and implement the veto which allows a 
satisfactory continuum limit at fixed Tq. Discussion and concluding remarks are the subject of the 
final Section. 



2 Discretized 't Hooft Model 



The Lagrange density for SU (Nc) gauge fields coupled to quarks in the fundamental representation 
is given by 

4 

where F^jy = d^A^^ — d^A^ — ig[A^, Ai,]. We remind the reader that the normalization of gauge 
fields appropriate for matrix fields and dictated by the gluon kinetic term differs by a factor 1 /\/2 
from the more standard one: 



C = --TtF^'T^, + q[i^.^d- igA) - i^o] q, (2.1) 



with Fg = ^-fa- Thus Ag = yl/\/2, and we conclude that g = gg/V^- In 2 space-time 
dimensions we choose the representation of 7 matrices for which the light-like components are 

With this choice the field equation for the upper component of the quark spinor does not involve 
the "time" derivative and is an equation of constraint relating the upper component, qi, to the 
lower component, q2- Working in light-cone gauge {A- = A~^ = 0), we can eliminate the upper 
component in favor of the lower component yielding the light-cone gauge Lagrange density 



£ = +^Tr{d_A+f + ii;^ 



i; (2.3) 



where ip = 2^/^(72- 

Our discretization of Feynman diagrams is based on the x"*" representation of each bare propa- 
gator 

D{p+,x+)= [ —D{p+,p-)e-'''^P~. (2.4) 
./ 27r 



Performing the p integral gives the following Feynman rules for the continuum theory 



DJp+,x+) = e-"-+'^o/2p+ ^ e-M^/V 



Da{p\.-) = ^^^-^ (2.5) 
p+ p+ 
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where the arrows indicate the rules to use with imaginary time. 



One way to digitize the 't Hooft equation (LI) is to put the variables x,y on a grid, which 
amounts to discrete light-cone quantization ^ , where one discretizes the amount of each line 
of the ladder diagram carries in quanta of m 



P 



Im 



I = 1,2,3,... . 



One can then focus on a state of the system of interest (in our case a qq system) with total 
P"*" = Mm. The continuum theory is recovered by taking the combined limits m — > and M ^ oo 
while keeping P+ = Mm fixed. Following ||l|, ^, in addition to discretizing the p'^ of each particle, 
we also discretize imaginary light-cone time, r = ix'^ = ka [k = 1,2,3,...). This discretization 
(which also serves as an ultraviolet cutoff) allows the continuum limit to be reached by keeping 
Tq = m/a fixed and taking both m, a — > and M oo simultaneously. Actually, since the physics 
of the discretized model depends only on the ratio m/a, the continuum limit is nothing but the 
large M limit, where M measures the total of the system state. The conventional continuous 
time DLCQ approach (see and references therein) corresponds to the special case Tq — > oo. 

Discretization of the quark propagator poses no difficulty. However, for the instantaneous 
interaction induced by integrating out A^, we allow for ambiguities as in |||]. The only constraint 



is that the discretized propagator become that of Eq. 2.5 in the continuum limit. This allows us 
to spread out the instantaneous interaction away from r = (see Q for further discussion). Thus 
the gauge propagator can be expressed as 



DA{Mm, -ika) = -fkj^ where 



k=l 



1. 



(2.6) 



We require that these arbitrary parameters fk rapidly vanish with increasing k. Using this dis- 
cretization, the Feynman rules for the discrete theory are summarized in Fig. |^. 
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Figure 1: Feynman rules for the discretized 't Hooft model. Discrete light-cone time flows up the 
page. 

For the purposes of this paper we shall not exploit the full generality of the set of {/fcj's. We 
restrict attention to the simplest version where the spread out interaction propagates only one unit 
in light-cone time, this corresponds to setting fi = l,/fc>i = 0. The Feynman rules of Fig. || 
can be further simplified if we absorb the negative sign from the anti-quark propagator into the 
corresponding vertex factor. We define new parameters 



a = e-'^o/2To and k^J^. (2.7) 

V 2vrro 

We also recall that in 't Hooft 's large Nc limit every additional pair of cubic vertices in the ladder 
sum corresponds to a completed color index loop, which produces a factor Nc- Thus we shall also 
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absorb a factor of \/Nc into each vertex. Simply put, all terms in the ladder sum are only dependent 
on the 't Hooft coupling g^N^- The simplified Feynman rules are presented in Fig. |2|. 
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Figure 2: Simplified discretized Feynman Rules for 't Hooft model. 



3 Single Time- Step Transfer Matrix 

Using the Feynman rules of Fig. ^ we can now proceed to set up a transfer matrix which evolves 
a singlet qq system one step forward in 2;"'"-time. Once the matrix has been determined as a 
function of the coupling, k, solving the eigenvalue problem will yield the bound state energies 
as functions of coupling. Since the scalar particle which mediates the Coulomb interaction only 
lives one time-step, any state can have at most two intermediate scalars. Thus for the simplest 
systems with jm = M = 3, 4, 5, 6 the number of states are 3, 7, 14, 25 (the number of states is 
(M - 1)(M2 - 2M + 6)/6 for general M). For illustrative purposes we shall explicitly present the 
transfer matrix for M = 4. 

For M = 4 there are 7 states namely: 

|3,1) =6^4|0), |2,2) = 44|0), |1,3) = &I4|0), 

|2,1,1) = 44410), |1,2,1) = 6|44l0), |1,1,2) = 6|44l0)> (3-1) 

ii, 1,1,1) =6l4'4io). 



where 6^, S ^ and are creation operators for the quark, anti-quark and intermediate gauge particle 
states (the subscript on these operators denotes jm). 

By construction each of the quark and anti-quark states has at least one unit oip^ jm. The 
matrix that evolves the system forward in can be factored into a matrix A that involves only 
propagators and a matrix V that involves vertices. Writing the state of the system as a column 



vector, T, with 7 components corresponding to the seven states in Eq. |3.l| , the transfer matrix 
equation is 

tT = AVT, (3.2) 

where 

A = diag ^Q^^/^ a, a^/^ -a^/s _ ^2/4^ _q,3/2^ ^2! ^ (3 3) 



4 



/ 1 








K 


K 








\ 





1 





— K 

















1 





— K 


— K 







K 


— K 


















K 





— K 



















K 


— K 















V 






















and the eigenvalue is t = e~ . Solving this eigenvalue problem will yield energy eigenvalues 
as a function of the coupling n. Note that the matrix AV is not hermitian, and because of the 
negative diagonal entries in A, the equivalent matrix \/^V\/^ is not hermitian either. Thus there 
will, in general be complex eigenvalues t. The best one can hope for is that the lowest lying energy 
eigenvalues (highest lying positive real part for t) are real. A satisfactory outcome for the continuum 
limit M ^ oo would be that the ground state energy and all the energy values with real parts of 
order 1/M above the ground state energy are real. Then the complex eigenvalues would be strict 
lattice artifacts. 

The existence of complex t eigenvalues is already evident at M = 4 as shown in Fig. ^, where we 
have chosen a = 0.5 which for definiteness we use in subsequent graphs unless otherwise indicated. 
The ground state (highest value) of t stays real and positive for all coupling. However the next 
2 excited states stay real only for k < Hc when they collide with eigenvalues that have emerged 
from t = (infinite energy) after which the eigenvalues become complex conjugate pairs. The hope 
is that for increasing M the number of lowest lying energy levels that remain real all the way to 
strong coupling should increase. For M = 4 analysis shows that the lowest energy eigenvalue (that 
of the ground state) stays well-separated from the other states (real and complex) for all couplings, 
see Fig. g We also see the eigenvalue solutions (again see Fig. ^) which are well behaved at weak 




Figure 3: Plot of the real solutions of i as a function of for the M — 4 single time-step transfer 
matrix. It is convenient to plot t rather than energy since then infinite energy corresponds to t = 0. 
Also note that the lowest lying states are those with the largest value of t. 
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coupling can merge with t = solutions (solutions which have t = at zero coupling correspond to 
infinite energy lattice artifacts) and become complex. Complex t solutions are not physical as they 
correspond to complex energies. This behavior is generic for our discretization, but as we shall see 
later, when the problem has been set up correctly, we can separate the lowest lying states which 
survive the continuum limit from the lattice artifacts. 




Figure 4: Plot of the lowest energy real solutions of E = — \n{t)/a as a function of for the M = 5 

and M = 6 single time-step transfer matrices. We see the appearance of additional real solutions 
at lower energies than the weak coupling ground state for k > Kc- 

However, when one performs a similar analysis for the M = 5 and M = 6 systems the lowest 
eigcnstate at weak coupling does not remain the ground state for all coupling. In both cases a 
complex solution at weaker coupling becomes real at larger coupling with a lower energy than the 
weak coupling ground state. Comparing this behavior for M = 5 and M = 6 suggests that for 
increasing AI this probably occurs at weaker coupling. Thus for large M the weak coupling ground 
state might only be valid for extremely weak (perhaps only infinitesimal) coupling. 

Conventional continuous time DLCQ corresponds in our discretization to — > since then 
To ^ oo. In order for our light-cone time discretization to be useful, the solution should work 
for all coupling (corresponding to all values of To). Here, in this single time-step analysis, we see 
that our most naive discretization does not satisfy this requirement. We shall have to modify the 
discretization in order to fix this. 

Since the continuum limit requires M — oo the single time-step analysis is also inefficient 
because the rank of the matrix to diagonalize is of 0{M^). However, as we shall show in the 
following section, writing the ladder equation in the form of a Bcthe-Salpeter equation (cxchangc- 
to-exchange rather than single time-step) will reduce the complexity of the eigenvalue problem to 
a matrix of rank of 0{M). 
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4 Bethe-Salpeter equation 



A more efficient way to solve the discretized 't Hooft model is by setting up a Bethe-Salpeter 
equation Q. Instead of a matrix equation that evolves the qq system one step forward in time, we 
can write down a system of equations (also a matrix equation) which evolves the system exchange to 
exchange. The simplification is that the intermediate state involves two (dressed) particles (M — 1 
possible states for general M) rather than two, three, and four bare particles as in the case of 
the single time-step transfer matrix. The trade-off is that the equations become more complicated 
because of the dressed propagators. 

In order to set up the Bethe-Salpeter equation it is necessary to work out the dressed quark 
propagator. In the context of this discretization the dressed quark propagator is just the sum of 
all possible iterated bubbles. There is no room for nested bubbles because fk>i = 0. While the 
bubbles extend only one time-step in x^, we must still allow for all possible P'^ routings through 
each bubble. 



Figure 5: Iterated bubbles which contribute to the quark propagator. 



The energy representation of the bare quark propagator carrying /m = I (without bubbles), 
obtained by multiplying by and summing over all /c > 0, is given by 



k=l 



ua 



i/i 



ua 



i/i 



1 — ua^/^ ' 



where u = 1/t = e . The contribution of a single bubble is 

i-i 



l/{l~r) 



r=l 



The full propagator is given by iterations of Eq. 4.1 and Eq. 4.2 as displayed in Fig. ^ 



k=0 



The denominator of the full propagator can be factored in two roots so that 



it-t+)it-t.)' 



where 



a 



i/i . 



1 ± Vl - Aa-^/^K^^i 



We can now partial fraction the full propagator 



(t+ - t. 



{t-t+) (t-i_) 



it, 



k=l 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



(4.6) 
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Expressing the full quark propagator as the sum in Eq. allows us to read off the time represen- 
tation of the full quark propagator for discrete r = ka. 

What we really need in order to set up the Bethe-Salpeter equation is a 'propagator' which 
propagates the qq system, including bubbles, between exchanges between the quark and anti-quark, 
see Fig. The 'propagator' which evolves the system forward between exchanges is then 




1 M-1 



Figure 6: Here is a section of the ladder sum between two exchanges. The quark on the left carries 
jm — I and the anti-quark on the right carries /m — M — I. 



M/l(M-l) °° / \ / \ 



k=l 



where s± are the roots for the anti-quark (obtained simply by replacing / in Eq. 4.5 by M — I). 
With some manipulation this can be simplified to 

VgM) = 5 ^ ^-^^ , 4.8 

where for brevity, we have defined 

= aV(A^-Os,. (4.9) 
We can now now set up the Bethe-Salpeter equations 

A/-/-2 2 M-l-1 2 



r=l r=l 
l-l 2 '-2 2 

^■^W = Y'^V^^{l-r)^.^{l-r) + Y,'^^::fl-r)^:^il-r). (4.10) 



r=\ r=l 



^ I, ^'i label two-particle states where the last ladder rung propagated forward in time from left 
to right or right to left, respectively. The first equation is graphically portrayed in Fig. |^. Since 
each of the quark/anti-quark propagators must carry a minimum of one unit oi jm there are 
only 2(M — 2) possible states 

^ \{V) : l</<M-2, (/) : 2 < / < M - 1. (4.11) 
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Figure 7: Parallelogram and trapezoid sections of the ladder sum. Internal variables label the 
number of units of jm carried by each leg. The quark and anti-quark propagators include 
self-energy corrections. 



Eq. 4.1C| is constructed by evolving the system from a state just after one exchange in the 
ladder sum to just after the next. The various P's in Eq. ^I.IO] correspond to the Feynman diagram 
contributions which are either parallelogram or trapezoidal sections which take a ^'h or ^'^i to a 
^'h or ^'ri. The parallelogram propagator sections are simply related to Eq. 4.5, 



Vgg{M-l)=Vgg{l). 



However, the trapezoidal segments must be independently determined 

PP(0 = Vgg{M-l), 



(4.12) 



(4.13) 



where 



(t+-t_)(s_ 



(1 - u^K^Y.\Yl^ 



M-U 



M-U 



1 



(4.14) 



In order to solve the matrix equation in Eq. |4.10| we would like to write it in the form of an 
eigenvalue problem yielding t = l/u as a function of n^. This is slightly complicated since the 
propagator segments involve S"s which appear together with factors of k^. By setting x = ''^^^ we 
can manipulate the equation to isolate t as the eigenvalue, with solutions tn{x)- This is achieved 
by rescaling ^'h, ^'^] by the denominator factor common to all P's, yielding 



M-l-l 



ai (1 - (1 - x^'m-i) ^' , (0 + E i^i^ra^'^^''''^^ " X^'m-I-.) ^'^ + 
(1 - x^S^S', 



t 



r=l 

M-l-2 



E 

r=l 



X 



'1 
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M-l-1 2 



r=l 



l-l 



ai (1 - x^'i) (1 - X^'m-i) ^' , (0 + E ^»i-r»'^^''''^ (1 - I - 

r=l 



r=l 



r=l 



(4.15) 



where a, = a^/'^^"^). 

This discretized equation has roughly twice the complexity of a straightforward discretization of 
the 't Hooft equation. The reason is that a rung propagating forward from left to right can couple 
to subsequent evolutions forbidden to a rung from right to left (and vice versa). See Fig. |8|, for the 
graphs responsible for this asymmetry. This is the reason we had to introduce a two-component 



Figure 8: Asymmetry in the densest configuration of exchanges in the same sense and opposite 
sense. The double arrow points to the effect of implementing the veto. 

Bethe-Salpeter wave function. An immediate consequence is that at k = each energy value is at 
least doubly degenerate, including the ground state. This feature is evident in Fig. ^ where the 
solutions of the BS equation are displayed for M = 4,5. All of the solutions seen in Fig. |3| are 
present, but in addition there are extra spurious solutions. For example, with M = 4, there is a 
second curve emerging from the k = ground state eigenvalue. For k > this extra eigenvalue 
curve lies below (in t) and well separated from the true ground level curve for all coupling. Similarly, 
for other values of M the Bethe-Salpeter method consistently reproduces all the solutions of the 
transfer matrix method, but it also adds spurious solutions due to the two-component nature of 
the wave function. 

One way to avoid these unwanted solutions is to slightly modify the discretized Feynman rules so 
that the rung will attach to the same lines whichever way the exchanged gluon propagates. As seen 
in Fig. ^, the asymmetry stems from the possibility of consecutive gluon emissions (absorptions) on 
immediately successive time steps. If this possibility is disallowed, the basic exchange rung can be 



taken to be the sum of the two different exchanges as in Fig. 10. In addition to removing unwanted 



solutions this veto rule also leads to simpler equations, with a more transparent continuum limit. 
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Figure 9: Real eigenvalues using the Bethe-Salpeter method for M = 4, 5. All solutions of the 
single time step method, see Fig. ^, are reproduced, but additional spurious solutions are present. 

As we shall see in the next section, it also produces a more physical strong coupling behavior than 
our original discretization. 





+ 




Figure 10: With successive emissions and absorptions vetoed, the two types of exchanges can be 
combined in a single rung. 



5 Bethe-Salpeter with Veto 

The Bethe-Salpeter equation for the discretized 't Hooft model, with the veto imposed as described 
at the end of the previous section, is 



*(0 



M-l-1 2 



r=l 



,.2 



(5.1) 



r=l 



where Vgg is defined in Eq. n^. After re-indexing both sums the equation can be written as 



^(/) = ^ ^^«iA+i/{M-.)p^_(,)^(,) + g _:i!!!_«i/r+i/(M-Op^_(,)^(,). (5.2) 
(£ r) 



M-l 2 



r=l+l 
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By imposing the veto we have reduced the rank of the eigenvalue problem from 2(M — 2) to M — 1. 

The new discretized equation is much easier to analyze in the formal continuum limit M ^ oo 
than the original. First define <I>(r) = 2?gq(r)^'(r), and rearrange Eq. |5.2| to read 



(5.3) 



M-l 

E 

r=l+l 



l-l 



(/ - ry 



l/l+l/(M-r) 



l/r+l/(Af-0($(^)_$(^))_ 



To formally examine the continuum limit we suppose that each discrete p'^ variable is large putting 
each I xM,| and take M ^ oo at fixed x. Then the right hand side of Eq. is set up to go to 
1/M times the r.h.s. of the continuum 't Hooft equation: 



A/-1 2 
r=l+l ^ ^' 



M 



-p 



dy 



l-l 2 

^{y) - $(x) 

{y - 



(5.4) 



Clearly, u must be chosen so that the l.h.s. is also of order 1/M. Next, it is easy to verify 
that = a^l^^^l^^'^^^^iji^ 1^ — \/l + Oilnl/t^)), so that the inverse propagator can be simplified, 
neglecting terms of order In M/M^, 



UlVqq{l) 



TT 



M 



1 1 - uiK^{iT^/6 - M/2l{M - I)) 



where we have defined ui 
now be simplified to 



6 2l{M-l)J uil + uiK^{iT^/6- M/2l{M -I))' 
ua^/^'^^~^\ The factor K multiplying <I> on the l.h.s. of Eq. 5.3 



(5.5) 



can 



a 



M/l{M-l) 



1 



^ a 



M/l{M-l) 



4 4 
K IT 

36 



4 4 
K IT 

36 



1 1 - uik^{tt^/6 - M/2l{M - I)) K^vr^ 
n't 1 + UiK2(7r2/6 - M/2l{M - I)) 

M 



+- 



6 l{M-l) 
1 1 — UiK^vr^/e k?tt'^ 

Ui \ + UlK^TT^/G 3 

M 



M 



1{M - I) 



+ 



l{M-l) 



K^TT^ 



6 



{l + UiK?TT'^/Qf 



(5.6) 



Now write u = uoe"^, where aA will be determined to be of order 1/M, so that ui = no(l + aA + 
{M/l{M - I)) In a) to order 1/M. Then no must satisfy 



/(no) 



4 4 
K TT 



1 1 — uqk^tt'^ k^tt"^ 



Un 



36 uqI + uokP-jt"^ 



0. 



(5.7) 



''Of course even for M large the equation does contain terms where I and M — I axe smaU (i.e. close to 1). In 
order for these contributions to not affect the solution to the continuum Bethe-Salpeter equation, the wavefunction 
must vanish at the endpoints. We shall see how this occurs when we evaluate the numerics later. 
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Then the continuum hmit reads 



1 



Mx 1 



In a + 



uof'iuo) 



(1 + U0K27r2/6)2 



p 



Muof'iuo) Jo 



dy 



(5.8) 
$(y) - $(x) 



The energy of the system is E = (lntio)/a + A, but the divergent first term is simply a physically 
irrelevant M independent constant, so it is consistent to identify P~ = A. Then Ai'^ = 2P^P~ = 
2MmA = 2MToaA. We also identify 



2ro<ilna + 



1 



and we obtain the continuum 't Hooft equation 

1 



X 1 



P 



uofim) Jo 

2-KUof'{uo) Jo 



(1 + noK^TrVe)^ 



dy 



<^iy)-<^{x) 

{y-xf 



(5.9) 



Comparing with Eq. 1.1, we see that the only effect on the continuum limit of keeping Tq finite is a 
finite renormalization of the gauge coupling — > —g'^/uof'{uo), and a coupling constant dependent 
shift in fi'^. Thus, the only requirement for identical continuum physics is that uof'{uo) be negative. 
Since a is a free parameter, we can access all positive values of fi'^ by tuning it. 

Eq. |5]^ implicitly relates uq to k via a cubic equation. Instead of solving this equation, it is 
more illuminating to use it to relate uq to the combination r] = uok^vt^/G 



no 



(l-??2)(i + ^) 



6ry 



6ry(l + r/ + 27]"^ 



We can also obtain the charge renormalization factor uof'{uo) in terms of rj: 



Uof'iuo) 



;i + r/ + 27/2)(l + ry + 77/2 



:i-r?2)2(i + ^)2 



the effective coupling in the 't Hooft equation 



2k^To 12r/(l-ry2)(l+77)To 



71" Uof'iuo) 

and the renormalized mass parameter 



7r2(l + f] -\- 7r/2 — rj^) ' 



12r?2(3 + r?2)To 

7r2(l + 7/ + 77/2 _ j^3^ 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



where we have used a = e~'^o/^"^°. 

As a check, note that the continuous time limit corresponds to Tq ^ oo or k? — > 0, whence 
uo — > 1 and 7] 0. Then the effective coupling Eq. ^.12 goes to 12ror//7r^ = 2Tok^ = g'^Nc/ir = 
glNc/2'K as it should. Next, with discrete time, we see that, in order to have real energy and k, 
(no > and 0), we must place the restriction < 77 < 1. Small k corresponds to small rj, 
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and large k corresponds to rj near unity. Interestingly, we note that the effective coupling in the 
't Hooft equation is small in both the small and large k regimes. 

It is easy to understand the small effective coupling at large k in terms of our discrete time 
Feynman diagrams. With discrete time, ^ oo causes the diagrams with a maximal number of 
powers of k? per time step to dominate. For example the qq propagator Dgg behaves in this limit 
as 

..^M/l{M-l) oo 
- , , ,y,y> = Y.(--')"'(^1^'m-i)'^ (5-14) 

SO that the propagator for 2k + 1 time steps is a^^/^^'^^~^\K^T,'{l]'^j_iY (K^7r^/6)^'^ in the con- 
tinuum limit. We see that away from the endpoints there is a factor of K^7r^/6 per time step in 
the continuum limit, which corresponds to each quark propagating exactly one time unit between 
interactions. Since this is the eigenvalue of the transfer matrix, we immediately infer the strong 
coupling value of uq = 6/K^7r^. Because of our veto, every exchange between quark lines occupies 
precisely two time steps and contributes only a single factor of k?. Thus each exchange costs a 
relative factor of in the strong coupling limit, and this relative factor is proportional to the 
effective coupling in the 't Hooft equation. More precisely, separating out the factor correspond- 
ing to the strong coupling propagation of the quark and anti-quark for two time steps, we have 
= (K^7r^/6)^(36/K^7r^), so the effective coupling for a single exchange is SG/k^tt^ for large k, in 
accord with the r/ — > 1 limit of Eqs. ^.10| , ^.12 . 



Now we turn to a numerical analysis of our discretized dynamics in order to understand how 
the continuum limit is approached in practice. As with the no-veto case in section |^ we can write 
this equation as an eigenvalue problem by rescaling ^ and isolating the eigenvalue i as a function 
of X = UK?. The resulting eigenvalue problem to solve is 



m) = 



1 - x'^Kf-i) 



+x 



^1 ^l/l+l/{M~r) 



r=l+l ^ 



(5.15) 



We use numerical procedures in MAPLE and MATLAB to find the eigenvalues t^ix) of the matrix 
on the right hand side of this equation as a function of x- The value of is different for each t„ 
since = xtn- However by varying < x ^ oo we can generate the real solutions, for all k^. 
In order to solve for complex i„'s we would need to vary x in the complex plane rather than just 
over positive real numbers. 

The problem of contamination of the lowest lying states by complex solutions has been solved 



by our veto prescription: The lowest lying state for M = 6 for Eq. 5.15 remains intact for all 



coupling K?, see Fig. 11, which should be compared against Fig. |^ where the lowest lying state was 



only the ground state for < 3. When we analyze Eq. 5.15 for increasing M (see Fig. 11 for 
M = 16, 32, 64) we see that the number of low lying states that remain uncrossed for all couplings 
increases with increasing M. We also see that the spacing between these states decreases as M 
increases. Recall that the solutions in Fig. |l^ have been generated for a = 0.5. 

In order to compare our numerical results for large values of M (hopefully close to the continuum 



limit) with the numerical results of 't Hooft we solve the Bethe-Salpeter equation in Eq. 5.15 
for = 0.5 and 

a = 1.16433 a = 1.04167 a = 0.70930. (5.16) 
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Figure 1 1 : Plots of the lowest lying energy eigenstates of the Bethe-Salpeter equation with the veto 
for M = 6, 16, 32, 64. Other states which occur at higher energies than those displayed have been 
omitted. 
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These three choices of a correspond to values of 't Hooft parameter, [i? = tx^ jg^^Nc-, taken to be 
0, 1 and 2.11^ respectively. These values of /i were used in [^. Fixing is equivalent to fixing r/, 
uq and gl^Nc/Ti, thus choosing a value for Jl determines a in Eq. 5.13 . 
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Figure 12: Plots of the three lowest lying states against 1/M. The three graphs correspond to 
choices of a and such that the continuum limit of = /g^^Nc = 0, 1, (2.11)^ so that we can 
compare these results against those of 't Hooft 

As we can see in Fig. [l^ plots of the three lowest lying energy levels against 1/M show curves 
that become linear with increasing M. These results can be fitted to the functional form 



aii; = ln(no) + |^exp(|| + ^), 



(5.17) 



where C2 and C3 parameterize the departure from 1/M behaviour away from large M. We used the 



data of Fig. 12 in the range 128 < M < 2048 to fit this equation. With the fitted value of ci we 
can calculate the mass square of the corresponding bound state. As discussed previously, the M 
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independent term in Eq. 5.17 is dropped in identifying P . Since 

^ 2P+P- = 2MTo{aE - Inuo), 

we have, for = 0.5, 

<2 _ 2ci 



0.22265 



+ ..., 



(5.18) 



(5.19) 



in units of q^^Nc/t:- The results of the fits are tabulated in Table |T] against the results of 't 
Hooft We see that for p, = 1 and 2.11, the results of our discretization match quite well those 





fl = 


't Hooft 


/1 = 1 


't Hooft 


/i = 2.11 


't Hooft 


ground state 


0.72 





7.25 


7.2 


24.23 


24.1 


1st 


7.57 


5.9 


17.26 


17.3 


38.17 


38.1 


2nd 


16.21 


14.3 


27.06 


27.2 


49.98 


49.8 



Table 1: Comparison of numerical fits for 128 < M < 2048 (for /i = we used 128 < M < 4096) 
in order to determine the bound state mass squared in units of g'^gNc/n for our discretized theory 
compared against the numerical results of the conventional continuous time approach of 't Hooft. 

of However, for /i = we increased the range of M to 4096, which still yielded a poor match. 
What we did note was that even for these sizable values of M, convergence for /i = is slow. When 
fitting the data for fi = for the ground state to Eq. 5.17 we are trying to force it to fit a coefficient 
to a 1/M term which is not supposed to be there. It is more appropriate to use the form 



a£; = ln(no) + ^exp + ^) , 



(5.20) 



where the power /3 of the leading behavior is fitted dynamically. We performed this refined fit 
to the three lowest lying states for /i = which yielded the results assembled in Table 0. These 





/3 


ground state 
1st 
2nd 


1.108 
1.025 
1.015 



Table 2: Table of fits to /? in Eq. 5.2C for /t = 



results provide numerical evidence that for /x = 0, the 1st and 2nd excited states do have a nonzero 
meson mass (i.e. the leading behavior is 1/M). However, the leading behavior for the ground state 
decreases more rapidly than 1/M and is consistent with zero meson mass. 

We next address the issue of slow convergence for /2 = by examining the form of the ground 



energy eigenvector for increasing values of M. It is well known that the solutions of Eq. 1.1 for 
H = do not vanish at the endpoints x = 0, 1; indeed the exact ground state is simply a constant. 



As we can see in Fig. 13, at finite large M the ground state solution of our discretized equation 
is ever smaller at the endpoints, and the progression of shapes is toward a more square profile. 
But even for M = 4096 the eigenvector has not yet converged to its limiting form. This should 
be compared with the solution for jl = 1 which rapidly approaches it's limiting form (see r.h.s. of 
Fig. fsl). We see that, for our discretized equation, the solution for the ground state decreases more 
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rapidly near the endpoints {x = and x = 1) as M increases, consistently with the shape eventually 
approaching a square profile at M — > cx). However, it is not hard to show that consistency of the 
continuum limit requires that the range in x over which the fall-off occurs must decrease less rapidly 
than \j\fM . This still allows an approach to a square profile but convergence is necessarily slower 
than one might have expected. In fact all solutions of the continuum 't Hooft equation with /i = 
have non-zero values at the end points. Thus we should expect slow convergence for all solutions of 
the /i = equation because the discrete solution tends to vanish at the endpoints but the limiting 
form does not. This effect does not occur for /i > because then the continuum solution vanishes 
at the endpoints, so a decent approximation to it can be achieved with relatively smaller M. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

X X 



Figure 13: Plot of the ground state eigenvector against x = l/M for increasing M for the cases 
/i = and 1. Each eigenvector is plotted for the range of M indicated in powers of 2. 



6 Discussion and Conclusion 

In this paper we have explored the efficacy of the discretization of large Nc QCD proposed in ^ by 
applying it to the well-understood 't Hooft model. For a smooth continuum limit over the whole 
range of bare coupling k, we had to introduce a refinement of the discrete time gluon emission 
vertex. This amounted to insisting that after an emission, at least 2 time steps had to intervene 
before the next emission, with a similar restriction on consecutive absorptions. In contrast, an 
absorption is allowed to immediately follow an emission and vice versa. With this refinement in 
place we found that the continuum 't Hooft equation describes the mass spectrum for all real k. 
However, the parameters that occur in the equation are renormalized from their bare values, as 
summarized in Eqs. 5.10| , 5.12, 5.13. 



An amusing outcome of this renormalization phenomenon is that the effective coupling goes to 
zero in both the small and large k limits. Perhaps this feature is a version of weak/strong coupling 
duality, much celebrated in recent developments in string/M theory. However, we must concede 
that 2 dimensional QCD may be too trivial to expect anything other than the usual continuum 
theory to emerge from any continuum limit. Another caveat against attributing much significance 
to this "duality" phenomenon, is that the physics of the continuum limit really only depends on 
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the ratio fi^ /NcQ^. This is because one can always choose the effective couphng as the fundamental 
unit of energy. Then the theories at different coupling but with the same value of this ratio (0 
for example) are physically identical: any differences in description can be removed by a change of 
units. 

At any rate, we conclude that the discretization of can be meaningfully applied to QCD 
in 2 space-time dimensions, with some intriguing hints about the nature of weak/strong coupling 
duality. An obvious and important limitation of the 2 dimensional case, however, is that the gluon 
has no dynamical degrees of freedom. Thus there is no opportunity for the P"*" of the system to 
be shared amongst an infinite number of gluons. This must occur for the fishnet diagrams to be 
relevant, and is allowed in higher dimensional space-time. The next step is to study the three 
dimensional case, the simplest gauge theory where fishnet diagrams can be relevant. 

Acknowledgements: We thank Klaus Bering for his helpful contributions in the early stages of 
this project. 
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